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When a linear chain of plasmonic nano-particles is exposed to a transverse DC magnetic field, 
the chain modes are elliptically polarized, in a single plane parallel to the chain axis; hence, a 
new chain mode of longitudinal plasmon-rotation is created. If, in addition, the chain geometry 
possesses longitudinal rotation, e.g. by using ellipsoidal particles that rotate in the same plane as 
the plasmon rotation, strong non-reciprocity is created. The structure possesses a new kind of 
chirality-the longitudinal chirality-and supports one-way guiding. Since all particles rotate in the 
same plane, the geometry is planar and can be fabricated by printing leaf-like patches on a single 
plane. Furthermore, the magnetic field is significantly weaker than in previously reported one-way 
guiding structures. These properties are examined for ideal (lossless) and for lossy chains. 

PACS numbers: 41.20. Jb,42. 70. Qs,78.67.Bf,42.82.Et, 71.45. Gm 

Keywords: plasmonic waveguide, sub-diffracting chain, one-way waveguide, longitudinal rotation, longitudi- 
nal chirality 



I. INTRODUCTION 

Linear chains of identical and equally spaced 
nano-particles were studied in a number of pub- 
lications [HQ- They support optical modes 
with relatively low attenuation and with no ra- 
diation to the free space if the inter-particle 
distance is smaller then the free space wave- 
length A, and then the modes can be much nar- 
rower than A. Hence the name "Sub-Diffraction 
Chains" (SDC). SDCs were proposed as waveg- 
uides, junctions, and couplers (J-Jl]. 

Recently, SDCs were suggested as candidates 
for one-way guiding Q. The physics there is 
based on creating an interplay of two types 
of rotations: geometrical chirality and Fara- 
day rotation (caused by longitudinal magne- 
tization). This interplay strongly enhances 
non-reciprocity and eventually leads to one-way 
guiding. The attractive features in Q are the 
3D nano-scale dimensions and the fact that the 
magnetic field bias is an order of magnitude 
weaker compared to some alternative structures 
Q. However, its chiral geometry is difficult to 
fabricate. Other works on one-way guiding em- 
ploy periodic structures of non-reciprocal ma- 
terial operating around their Braag point, thus 
requiring relatively large structures [9|. 

Here we suggest a new type of SDC for one- 
way optical guiding. The structure is planar 
in nature, so it can be fabricated by relatively 
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simple printing procedures of leaf-like patches 
on a single plane. Moreover, it requires magne- 
tization field that is significantly weaker than in 
@, H[ ■ The underlying physics and geometry are 
described in Figs. [T][2j A conventional SDC of 
spherical plasmonic particles supports three in- 
dependent electric dipole modes; one with lon- 
gitudinal polarization p z , and two degenerate 
modes with transverse and mutually orthogonal 
polarizations p x , p y [J]. If this SDC is exposed 
to transverse magnetization, as shown in Fig. [TJ 
the p v mode is unaffected, but Bq = "QBq 
couples p x and p z . This coupling creates two 
new modes of elliptical polarizations in the x, z 
plane, with two new dispersion curves (see anal- 
ysis below). At any operation frequency, one 
wave propagates in the direction of rotation 
("paddles" forward and "rides" forward) and 
the second propagates counter to the rotation 
(paddles forward and rides backward). Hence, 
these new SDC modes are non-reciprocal. How- 
ever, their dispersion is still reciprocal (i.e. even 
in j3). The route to enhanced non-reciprocity 
and one-way guiding is to create an interplay 
of two-type rotations. Hence, we add a longitu- 
dinal chirality: we replace the spheres by ellip- 
soids, rotated in the same plane of the elliptical 
polarization, i.e. in the x, z plane as shown in 
Fig. [5] The rotation step is A9. As we show be- 
low, this chain indeed supports one-way guid- 
ing. It is periodic only for rational A8/n. If 
A8/tt = n/m and m,n are coprime, then the 
period consists of m particles. Unlike the spi- 
ral structure in Q, if this ratio is irrational 
then there is no coordinate transformation un- 
der which the chain becomes periodic. 
Below we use the Discrete Dipole Approxi- 
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FIG. 1. Top: a chain of spherical plasmonic par- 
ticles with transverse magnetization. This chain 
supports modes with longitudinal rotation. Bot- 
tom: supported modes dispersion. 




FIG. 2. A chain of plasmonic prolate ellipsoids with 
transverse magnetization and longitudinal chirality. 
This chain supports one-way optical guiding. 



mation (DDA) and polarizability theory. These 
are standard tools in SDCs analysis [t |(f. 
They hold when the particle radius a and inter- 
particle distance d satisfy a <C A and a <C d, 
but studies show good accuracy even for d = 3a 
10]. 



II. FORMULATION 

If a small particle with polarizability a is 
subject to an electric field whose local value in 
the absence of the particle is E L , its response 
is described by the electric dipole p = aE L . 
The dynamic tensor-polarizability of a general 
ellipsoid whose principal axes are aligned with 
the x,y,z axes, made of an anisotropic material 
with electric susceptibility x^ is obtained via 
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Here k is the vacuum wavenumber, I is 
the 3x3 identity matrix, and the imagi- 
nary term i£; 3 (67r) -1 I represents radiation loss. 
V = 4:ira x a y a z /3 is the ellipsoid volume with 
a x ,a y ,a z its principal semi-axes, and L = 
diag(iV x , N y , N z ) is the depolarization matrix 



whose entries are obtained from a x , a v 



by el- 



liptic integrals and satisfy J2 U N u — 1 Un- 
der the Drude model and magnetization Bq = 
yB , X is 
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law, Xyy = Xxx ~ 

- iujuib, and where 
and where uj p and 



ia), and Xxz ; 

L0/uJp, Qb = u b /uj p , 

uji, = —q e Bo/m e are the plasma and cyclotron 
frequencies, a = (ra p ) _1 represents loss, where 
t is the dissipation time constant. Equations 
([1])- fully describe the ellipsoid at the origin. 
The n-th particle polarizability, oc n , is given by 
the nA9 longitudinal rotation (about y) 
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where the non-zero entries of T„ are t\\ = £33 = 
cosnA(9, £22 = 1, £13 = — *3i = sinnA#. The 
electric field at (0, 0, z) due to a short dipole 
p at (0, 0, z') is given by the matrix relation 
E{z) — 6q 1 A(z — z')p where 
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here Ai = diag(l,l,0), A 2 = diag(-l, -1, 2). 
We express now the local exciting field of the 
m-th particle in the chain as a sum of contri- 
butions from all its neighbor-dipoles, and apply 
OL m . The result relates the m-th dipole excita- 
tion p m to its neighbors 

Pm = e Q la m A t( m _ ^ Pn- ( 5 ) 

n, n^m 

This equation is not shift-invariant. Further- 
more, since is a rotation about y and not 
about the longitudinal axis z, it does not com- 
mute with the propagator A so the mathemat- 
ical transformation used in cannot be ap- 
plied. However, we set A9/n — N/M rational, 
hence the period D = Md consists of M parti- 
cles with polarizabilities c*o, £*i, ■ • ■ £*m-i as in 
Eq. Q and a m +iM — Q m V integer I. Then by 
periodicity 
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and the chain modes are determined from the 
M vectors of the reference period p , . . . , Pm-i- 
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We substitute this solution into Eq. ([5]). For 
each to within a period we decompose the in- 
finite sum into a sum Sq of contributions from 
identical particles (n = to + IM, 1^0), and a 
set of summations S m - n ' of contributions from 
the rest (n = n' + IM, n' ^ m). The result is 
the 3M x 3M matrix equation 

M-l 

(foa™ 1 - So)P m - J! S ™-nP n = 0, (7) 
n— O.ny^m 

with m = 0, . . . , M — 1, and where 

(E^ A(lD)e^ D , q = 0, 
S q ={ (8) 
[EiMqd - W)e^ lD , l<|g|<M 

The modes are obtained by looking for w(/3) 
at which Eq. (0 determinant vanishes. The 
vectors that span the corresponding nullspace 
are p , . . . , Pm-i that with Eq. ([6]) describe the 
entire chain excitation. 

Note that the series for S q above converge 
poorly. However, they can be cast in terms 
of the Poly logarithm functions Li Sl for which 
efficient summation formulas exist (see [1, [l3[ 
and Appendix). First observe that the matrix 
Sq is identical to that obtained in conventional 
chains with d = D. Hence [4[ 

P 3 

S = —J2uJs(kD,j3D)A s , (9) 

s=l 

where (ui, 112,113) = (1, —i, 1), A3 = A 2 and 
f s (x,y) = x- s [Li 8 (e ix+i y) + Li^e^)} (10) 

where Li s (z) = Y^=i 7^ ^ s the s-th order Poly- 
logarithm function. To write S q in terms of Li s , 
note that all sums in S q have the general form 

« ( e iC)lM+ g ~ / 
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where a„(q) = a n+ M{q) is M-periodic sequence 
with a n (q) = 0forl<n^g< M, and a q (q) 
1 . But M«„(g) = X^io 1 e 227rr ( n -«)/ M , hence 

Af-1 

C = Af" 1 X] e" i2,rr « /M Li a (e* +i2,rr/M ). (12) 

r=0 

Following the standard steps leading from 
Eq. © to Eq. © and using Eq. ([P2")). we get 
for S 9 , 1 < |g| < M the same expression as in 
Eq. but with f s replaced by h s (kd, j3d;q), 

(13) 



III. EXAMPLES 

We turn to study the structure. First, we ex- 
amine a transversely magnetized chain of spher- 
ical particles of radius a. It is a special case of 
Eq. ([7]) with the sum over S q ,q 5^ dropped, 
M = 1(D = d), A6 = 0, and a.Q = a given 
by Eqs. CO <(2]) with all depolarization factors 
set to 1/3. The chain dispersion is obtained 
by Det(eoa _1 — So) = 0. For lossless chain, 
this can be satisfied only if Im[diag(S'o)] = 
— fc 3 /(67r). Also, the first and third rows of 
Eq. now are linearly dependent, and read 

{w 2 -l/3 + y[(5 )ii+ifc 3 /(67r)]}p x = iujQ bPz . 

However, as discussed above the term in the 
square brackets must be real. Hence there is a 
phase difference of ±7r/2 between p z and p x if 
uii, ^ 0. This implies that the chain modes are 
elliptically polarized in the (x, z) plane. 

The chain dispersion for d — 3a = A p /30 and 
no loss was calculated numerically, and is shown 
in Fig. [TJ We have applied a relatively strong 
magnetization of uji, — 0.05a;p, in order to ob- 
serve the features clearly (the one-way property 
shown below is obtained at much weaker mag- 
netizations). The y-polarized mode is identical 
to that of a conventional chain. Two additional 
modes of elliptic polarization in the x — z plane 
exist in the chain. When observed from y > 0, 
one rotates clockwise [shown by the upper (red) 
curve] and the second rotates counter-clockwise 
[shown by the lower (black) curve] . This elliptic 
polarization is only due to the transverse mag- 
netization. It is mathematically evident from 
the ±7r/2 phase difference between the x and 
z components we have observed in the solu- 
tions for p , as seen in Eq. (fill) . All dispersion 
curves are even in j3, permitting propagation in 
both +z and —z directions. Also all dispersion 
curves possess the light-line gap clearly seen in 
the center. However, the light-line dispersion 
branches that run parallel to the light-line cone 
j3 = lo 1 1 c and that are associated with transverse 
polarization in conventional chains [U, [H| , rep- 
resent modes that are practically not excited 
[TB ] in the elliptically polarized curves (red and 
black). To avoid cluttering the figure they are 
not shown here. 

Next, we add longitudinal chirality; we re- 
place the spheres by prolate ellipsoids with 
semi-axis a x = A p /90, axes ratio of a y — a z = 
0.9a x (nearly spheres), and add longitudinal 
chirality with A9 = 7r/3 (hence chain period 
is D = 3d = Ap/10). Figure [3] shows the 
results for the relatively weak magnetization 
of lui, — 0.005oj p . Due to the added chiral- 
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ity, the upper (red) dispersion curve of Fig. [T] 
that corresponds to elliptical polarization with 
clockwise longitudinal rotation, splits to three 
branches shown in Fig. [3^,, and shifts the light- 
line gap rightward. Hence symmetry is broken 
and one way guiding is supported at frequen- 
cies within the shifted gap, as shown in Fig.|3Vj. 
Finally, in Fig. [3J; we show the normalized re- 
sponse \p n \ I |pq| of a finite chain of N = 600 
particles with the parameters above, for a unit 
amplitude excitation of the central particle (at 
the origin), at the selected frequency. This 
response is obtained by solving Eq. ([5]) for a 
finite number of particles, with only the cen- 
tral particle forced to oscillate at the aforemen- 
tioned frequency (hence it is a matrix equation 
of 3(N — 1) unknowns). The three curves in 
Fig. [3J: correspond to the three particles within 
each period. One way guiding is clearly ob- 
served (note the logarithmic scale) . Essentially 
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FIG. 3. One-way guiding in the upper band, (a) 
Dispersion curves. The gap in Fig[T]is shifted right- 
ward by the longitudinal chirality. (b) Frequency 
selection for one-way guiding, (c) Chain response. 

the same picture holds for the lower (black) dis- 
persion curve of Fig. [I] that corresponds to el- 
liptical polarization with counter-clockwise lon- 
gitudinal rotation, when the above longitudi- 
nal chirality is added to the structure. This is 
shown in Fig. |H 





We note that these results are obtained with 
a very slight breach of spherical symmetry 
(a z = 0.9a x ). This low ellipticity is sufficient to 
create one-way guiding if Bq = yB^ is present. 
In addition, Bq is considerably weaker than 
that required in previous studies on one-way 
plasmonic waveguides. 

Next, we study chains with loss. Since a real 
dispersion ui(/3) does not exist, we have solved 
Eq. ((SJ) for a finite chain with N = 250 and 
for the same geometrical parameters (AO, d) as 
in the lossless example of Fig. 01 The operation 
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FIG. 4. One-way guiding in the lower band. 



FIG. 5. One-way guiding in a lossy chain, (a) 
Chain response, (b) Isolation ratio. 

frequency w and axes ratio r = a z /a x were cho- 
sen by scanning over the neighborhoods of uj, r 
for which the corresponding lossless chain pos- 
sesses one-way property. The results are shown 
in Fig. [5l Figure shows the chain response 
for loss parameter a = (rw,,)" 1 = 8.2 x 10~ 4 
corresponding to Palladium (Pd) [3] , and mag- 
netization of u>b = 7.5 x I0~ 3 ujp. It is obtained 
with w ~ 0.57cl>p and r = 0.92. For a clear as- 
sessment of the one-way property as a function 
of magnetization and loss, we plot in Fig. \Ejp 
the "isolation ratio" |p_„|/|p„| vs. n for var- 
ious values of loss and magnetization. Again, 
one-way guiding is evident. 

In the examples above the particles are nearly 
spheres. These examples are important from 
the theoretical/physical point of view. They 
show that by using our approach of two-type 
rotation interplay a very slight breach of sym- 
metry is sufficient to create a profound non- 
reciprocity and one-way guiding. However from 
the practical point of view structure flatness is 
at least as important, since it opens the way 
to employ planar fabrication technologies. To 
demonstrate this option we have simulated a 
chain with the following properties. Axes ra- 
tio is a x : a y : a z — 1 : 0.1 : 0.75, hence the 
particle's size in the y-direction is an order of 
magnitude smaller then its dimensions in the 
x, z plane. Since the longitudinal chirality is 
obtained by rotating the particles about their y- 
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axis, it follows that the entire structure is prac- 
tically flat and lies in the x, z plane. Other geo- 
metrical parameters are as before, with magne- 
tization level of LJb = 0.0075w p . A strong one- 
way property exists and it is shown in Fig. [6] 
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FIG. 6. The same as Fig[3] but for flat geometry 
with a x : a y : a z — 1 : 0.1 : 0.75. 

Finally, we verify our approach via a full- 
wave numerical solutions using the CST soft- 
ware package for the flat structure shown 
in Fig. [7J The particles are disks of 2.5nm 



apply stronger magnetization. Hence we use 
ujb = 0.03uj p . We have used the DDA for a first 
estimate of the chain operation parameters, and 
then used the CST to search more accurately 
around the initial guess and to get full- wave so- 
lutions. The full- wave based chain response as 
a function of frequency is shown in Fig. [SJ One- 
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FIG. 8. E-fields just above the upright particle in 
periods#±10 and ±15, obtained by full-wave so- 
lution of the chain shown in Fig. [7] For each fre- 
quency, the values are normalized to the field above 
the upright particle adjacent to the dipole from left, 
shown in inset (b) of Fig. [7] 
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FIG. 7. A flat chain geometry for one-way guiding, 
(a) Geometry details, (b) Excitation by a short 
dipole located in between two identical chains. 

thickness and of elliptical shape in the x, z 
plane with long and short diameters of 25nm 
and 20nm, respectively. The inter-particle dis- 
tance is d — 27.5nm. The longitudinal chirality 
angle is A8 = —60°. The particles material 
plasma frequency ui p and loss parameter t _1 
are 2tt x 2000THz and 2THz, respectively (cor- 
responding to Ag). Due to a limited availabil- 
ity of computing power and memory resources, 
we can apply the full-wave solution to chains 
of about 100 particles or less. Hence our chain 
consists of 92 particles that create two identical 
chains of 46 particles each, one chain on each 
side of a short dipole antenna that serves as 
a local source - see inset (b) in Fig. [7J Since 
this geometry should show significant one-way 
guiding within a distance of 46 particles only 
(instead of within hundreds of particles used 
to obtain the results of Fig. [SJ we need to 



way behavior can be clearly observed within the 
frequency band of 541-543 THz. The ratio be- 
tween the chain excitation on left and right side 
of the chain-the isolation ratio-is in the order 
of 20dB. Finally, Fig. U shows the E-ficld in the 
y = plane at the central frequency of the one- 
way band, 542THz. 




FIG. 9. Full-wave solution for the E-field in the 
y — plane, at 542THz. Normalized fields in the 
center of each upright particle are shown on top. 

Note that in this frequency A ~ 0.55/im. 
Hence a significant one way behavior is clearly 
observed over distances of 0(A). 



IV. CONCLUSION 

When an interplay of two-type rotations is 
supported in a guiding structure, strong non- 
reciprocity and one-way guiding are created 
Q. In this work, a new geometry possessing 
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longitudinal chirality in sub-diffraction parti- 
cle chains is proposed. When combined with 
transverse magnetization, the resulting ellip- 
tical rotation of the chain dipole-modes and 
the geometrical rotation (longitudinal chiral- 
ity) provide the two-type rotation that lead 
to strong non-reciprocity and one-way guid- 
ing. The required magnetization is consider- 
ably weaker than in previously reported stud- 
ies. This one-way chain can be fabricated by 
thin-layer printing of leaf- like flakes of thin met- 
als. Since single-atom layers of Graphene may 
behave as thin metal flakes with controllable 
properties, exhibiting SPP like behavior with 
very low loss [l6| , the present work may offer a 
basic scheme for one-way guiding on Graphene 
layer. This is a subject of ongoing research. 



Collecting the Clausen integrals components 
that constitute our Li s and rearranging, we ob- 
tain for s = 2 and Ixl < 7r 



Li 2 (e ix ) 



7T 



X 



6 4 

x In x — x 



(2tt - x) (A2) 



IE 
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where £?2rt are the Bernoulli numbers absolute 
value. With terms up to n = 5 in the series 
above, the relative error for x = it is in the order 
of 10 -5 . From the integral relation Eq. (IA1[) we 
obtain 



Appendix A: Evaluation of Li 3 (e t6 ) 

Li s {e lx ) for real x, can be expressed as fi- 
nite sums of Clausen integrals . p 7 . The 
latter can be computed using series that con- 
verge much faster than the algebraic conver- 
gence of the formal definition in Eq. ([5]) . Note 
that Lii(z) = — ln(l — z). Furthermore, from 
Eq. © it follows that 

Li s+1 (z)= [ r 1 Li s (t)dt (Al) 
Jo 



Li 3 {e™) =(,{?,) + —(\nx-Z/2) 

4 ^ n(n + l)(2n + l)(2n)\ 

71=1 

. / 7T 2 X 7TX 2 X 3 \ , . , 

+ \— ~ — + T2) (A3) 

where C(3) « 1.2020569. Higher order Polylog- 
arithm functions are not needed in the present 
work. 
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